OPEN 3 ACCESS Freely available online 

Effects of Caste on the Expression of Genes Associated (§S 
with Septic Injury and Xenobiotic Exposure in the cros^rk 
Formosan Subterranean Termite 

Claudia Husseneder*^, Dawn M. Simms^ 

Department of Entomology, Louisiana State University Agricultural Center, Baton Rouge, Louisiana, United States of America 

Abstract 

As social insects, termites live in densely populated colonies with specialized castes under conditions conducive to microbial 
growth and transmission. Furthermore, termites are exposed to xenobiotics in soil and their lignocellulose diet. Therefore, 
termites are valuable models for studying gene expression involved in response to septic injury, immunity and 
detoxification in relation to caste membership. In this study, workers and soldiers of the Formosan subterranean termite, 
Coptotermes formosanus, were challenged by bacterial injection or by no-choice feeding with a sublethal concentration 
(0.5%) of phenobarbital. Constitutive and induced expression of six putative immune response genes (two encoding for 
lectin-like proteins, one for a ficolin-precursor, one for the Down syndrome cell adhesion molecule, one for a chitin binding 
protein, and one for the gram-negative binding protein 2) and four putative detoxification genes (two encoding for 
cytochrome P450s, one for glutathione S-transferase, and one for the multi antimicrobial extrusion protein), were measured 
via quantitative real time polymerase chain reaction and compared within and among 1) colonies, 2) treatment types and 3) 
castes via ANOVA. Eight genes were inducible by septic injury, feeding with phenobarbital or both. Colony origin had no 
effect on inducibility or differential gene expression. However, treatment type showed significant effects on the expression 
of the eight inducible genes. Caste effects on expression levels were significant in five of the eight inducible genes with 
constitutive and induced expression of most target genes being higher in workers than in soldiers. 
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Introduction 

As social insects, termites (order Blattaria, formerly Isoptera [1]) 
live in densely populated colonies with a specialized caste system 
characterized by morphological differentiation and division of 
labor between worker, soldier and reproductive castes. Related- 
ness within termite colonies is usually high since members of all 
castes are the offspring of a limited number of reproductives (kings 
and queens [2]). Morphological and behavioral differentiations 
among castes are therefore based on gene expression rather than 
genotype [3]. 

Termites must have evolved strong and versatile defense 
mechanisms against toxic xenobiotics and pathogens in adaptation 
to living in densely populated colonies of closely related 
individuals, with intense social interactions in nest conditions 
conducive to microbial growth and transmission as well as being 
exposed to contaminants in the soil, insecticides and a lignocel- 
lulosic diet containing secondary metabolites produced by plants. 
Detoxification enzymes to metabolize xenobiotics are therefore 
essential. Previous studies using purification and biochemical 
characterization showed that termites possess multiple cytochrome 



P450s [4], glutathione S-transferases (GSTs [4]) and microsomal 
esterases [5] and show different levels of susceptibility to 
termiticides [6,7,8]. However, the relationship between constitu- 
tive or inducible gene expression and enzyme activity/insecticide 
susceptibility has not yet been established. 

While detoxification enzymes protect insects from toxic effects 
of xenobiotics, the innate immune system is the first line of defense 
against microbial infections. The innate immune system is 
primarily non-specific and serves to recognize invading pathogens, 
activate humoral and cellular immune responses leading to 
inactivation and clearance of pathogens. 

The insect innate immune system recognizes invading patho- 
gens through an array of pattern recognition receptors (PRR) that 
bind to conserved microbial cell membrane or wall epitopes 
serving as immune elicitors, such as lipopolysaccharides (Gram- 
negative bacteria), peptidoglycans (Gram-positive bacteria) and P- 
1,3-glucan (fungi). Insect PRRs include C-type lectins, gram- 
negative bacteria binding proteins, peptidoglycan-, lipopolyssac- 
charide-and P-l,3-glucan binding protein, and P-l,3-glucan 
recognition protein [9,10,11]. Binding of microbial immune 
elicitors by PRRs triggers cellular and humoral immune responses. 
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Cellular defense consists of a variety of hemocytes responsible for 
sealing wounds via aggregation [12], phagocytosis, encapsulation 
of invading microbes and nodulation, which leads to melanization 
by the prophenol oxidase (PPO) pathway and sequestration 
[13,14]. Humoral response includes the release of antimicrobial 
peptides (AMPs) into the hemolymph mainly from the fat body 
and the activation of proteases, which trigger melanization via the 
PPO pathway ([15], and references therein). The PPO pathway 
also produces intermediates that are direcdy toxic to microorgan- 
isms (e.g., quinones). 

While the PPO pathway initiates melanization and pr<)duc(;s 
antimicrobial compounds within minutes, recognition of foreign 
cell wall components also triggers intracellular signahng pathways 
that work within a couple of hours to days. The Toll pathway 
responds mainly to Gram-positive bacteria and fungi, while the 
immune deficiency (IMD) pathway reacts predominantly upon 
infection with Gram-negative bacteria. Both pathways result in 
expression of PRRs and antifungal peptides and/or AMPs. The 
Toll pathway also triggers activation of serine proteases involved in 
melanization. Other signaling pathways [P»IK- and JAK/STAT) 
are triggered by exposure to lipopolysaccharides of Gram-negative 
bacteria, but their contribution to insect immunity is less well 
researched [16]. 

Termites possess both humoral and cellular immunity (reviewed 
in [17]). Pioneering work by, e.g., [18,19] showed that immuni- 
zation via sub-lethal dosages of pathogens induces the production 
of AMPs in the termite hemolymph. A number of genes and 
peptides related to this kind of humoral immunity ha\'(; been 
described in termites. Among them are a transferrin gene which 
was up-regulated following exposure to an entomopathogenic 
fungus [20], termicin, spinigerin [18] and gram-negaUve binding 
proteins (GNBPs) which can act as both PRRs and AMPs in 
several termite species [21]. Production of AMPs appears to differ 
not only between species [22], but also between termites belonging 
to different colonies or castes [17,19,23]. The components of 
cellular immunity in termites are less well studied although 
melanization and encapsulation of foreign bodies have been 
frequently used as bioassay for immune competence (e.g., [24]). 

While both, detoxification enzymes as well as immune proteins 
have been detected in termites, full genome screening and gene 
expression studies are still in their infancy where termites are 
concerned. We chose the Formosan subterranean termite (FST), 
Coptotermes formosanus Shiraki (Blattodea: Rhinotermitidae), as 
model to investigate constitutive and induced gene expression 
related to xenobiotic metabolism and immunity in termites. This 
termite represents one of the most successful and economically 
important invasive pest species worldwide [25,26]. As a subterra- 
nean termite, C. formosanus forages and nests mainly in the soil 
and is associated with urban environments and as such exposed to 
xenobiotics not only through their natural diet (wood) but also 
through soil cxmtaminants and insecticides. Colonies are densely 
populated and can contain over a million individuals [2 7] . Such a 
hfe style should be conducive for growth and spread of pathogens, 
however, 50 years of failed biocontrol efforts [28] suggests the 
existence of strong synergistic immune and detoxification defense 
mechanisms of subterranean termite colonies. 

Annotation of our previously constructed expressed sequence 
tag library of the FST revealed the presence of putative 
detoxification and immune-related genes [29]. Out of those, we 
selected the 10 genes with the most similarity to known 
detoxification genes (e.g., those encoding cytochrome P450s, 
glutathione S-transferase) and immune genes (e.g., those encoding 
for gram-negative binding proteins, lectin- and fibrinogen-like 
peptides as well as the Down syndrome cell adhesion molecule 



[29]). A recent study aiming at a more complete description of the 
immune related transcriptome in C. formosanus [30] confirmed 
the presence of those genes among hundreds of putatively 
immune-related genes encoding for PRRs, signal modulators, 
transductors and effector molecules (such as AMPs) whose 
expression was upregulated after fungal and/or bacterial chal- 
lenge. However, colony or caste specific gene expression, inducible 
detoxification or possible crosslinks between genes involved in 
immunity and xenobiotic metabohsm were not considered. In fact, 
very few studies, have considered differential gen(; expression of 
immunity and detoxification genes among cast(\s of tc'rmites 
[3 1,32] . Therefore, the principal goals of our study were to identify 
genes associated with immunity and/ or detoxification in FSTs, test 
whether expression can be induced by septic injury with bacteria 
and/or sublethal xenobiotic challenge with phenobarbital and to 
test if expression of immune response and/ or xenobiotic metab- 
olism genes is colony or caste dependent. 

Results 

1. Constitutive and induced gene expression 

Examined in this study were a total of 10 genes putatively 
involved in immunity, and/or detoxification in termites (Table 1). 
Constitutive (untreated controls) expression of target gene 
transcripts was rather low, overall (mean R= 0.014, SD = 0.163, 
Table SI). AU but two target genes were inducible, either with a 
sublethal dose (0.5%) of the xenobiotic phenobarbital (PB) (mean 
R = 2.014, SD= 1.379) or septic injury with non-pathogenic 
bacteria {Escherichia coli/Pilihacler lermitis) (mean R = 2.668, 
SD= 1.988) or both treatments (Table SI, Figure 1). Only one 
of the cytochrome P450 genes (CYP450) and the gene for Down 
syndrome cell adhesion molecule (Dscam) were not induced by 
either treatment. Most genes were inducible by both treatments. 
The genes inducible by xenobiotic and septic injury/immune 
challenge encoded for CYP15A1, Chitin-binding protein (CBP), 
two lectin-related proteins (LECT-like, CLECT), the Ficolin-2 
precursor (FICO-2) and the gram-negative binding protein 2 
(GNBP-2). The gene encoding for glutathione S-transferase (GST) 
was exclusively induced by 0.5% PB treatment, while the gene for 
the Multi antimicrobial extrusion protein (MatE) was induced 
predominantly with the E. coli/P. termitis treatment. 

2. Effects of colony origin 

Tests for variance between colonies and for interactions 
between colony and treatment or between colony and caste 
showed no significant effect of colony origin at any of the 1 0 target 
genes. Probabilities of colony variability and interaction ranked 
from 0.127 to 0.783. Thus, regardless of treatment or caste, R 
values were not significantiy different among the three termite 
colonies. 

3. Effects of treatment 

In contrast to colony origin, treatment type (0.5% PB vs. E. coli/ 
P. termitis) had a significant effect on the expression of each of the 
eight inducible target genes (Figure 1, Table 81). CYP15A1 was 
induced in both types of treatment; however, relative expression 
was significantly higher (F = 88.03, P<0.0001: PROC GLM) for 
the E. coli/P. lermitis treatment as compared to the 0.5% PB 
treatment due to a strong upregulation of gene expression in 
workers after bacteria treatment. Similarly, GNBP-2 showed 
induced expression in both treatment types but with significantiy 
higher expression (F = 28.08, P<0.0001) for termite workers 
injected with E. coli/P. termitis than for termites fed with 0.5% 
PB. Also, CBP was significantiy more induced (F= 31.19, P< 
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Table 1. Descriptions of target and reference genes. 





Gene 


Abbrev. 


Species 
of Origin 


Top Matches in GenBank, 
Species of Origin 


Putative 
Function 


Primer 

Sequence (5 -3 ) 






(Accession No.) 


(Accession No. for Nucleotide 
or Peptide Sequence) 






Putative xenobiotic metabolism related genes 


Cytochrome 
P450 15A1 


CYP15A1 


Coptotermes 
formosanus 


CYP151A, Reticulitermes 
/7aWpes(ACN93795) 


Epoxidase, 


ATTTGTCACTCGCCTTGGTC 






FK835449 


CYP151A, Diploptera punctata 
(AAS 13464} 


juvenile hormone 
synthesis 


AAAGTTGCCATACGTGGAGG 


Cytochrome 
P450 


CYP450 


Coptotermes 
formosanus 


CYP314a1, Drosophila 
melanogaster (Q9VUF8.3) 


Ecdysone 20- 
monooxygenase 


ACTGGTTTGTAGTCGGTGCC 






FK833823 


CYP12A2, Musca domestica 
(018635.1) 


Monooxygenase, 
oxidoreductase 


TTAGCACCGAGAGACGTTCA 


Glutathion-S- 
Transf erase 


GST 


Hodotermopsis 
sjoestedti 


Glutathione S transferase, 
Blattella germanica 


Metabolism of 
lipophilic toxins, 


GTTTCAAGCCTCGCTTTTGAC 






DC239424 


(FJ855500.1) 


excretion of 
xenobiotics 


CGTCACTGTGAGACTGCCAT 


Multi 

antimicrobial 
extrusion 


MatE 


Coptotermes 
formosanus 


MatE, Trichonomonas vaginalis 
(XP_001314814.1) 


Sodium 
antiporter, 
cationic efflux 


TTGACACAACACCTTTCGGA 


protein 




FK833694 




transport pump, 
excretion of 
organic cations 


GCCGTGTTTGAGGATGAAAT 


Putative 
Immune 
related genes 


Chitin binding 
protein 


CBP 


Coptotermes 
formosanus 


chitin-binding protein, 
Holotrichia oblita 
{AFD28282.1) 


Peritrophin A-type 

chitin-binding 

protein 


CGGGAGGACTGCTAGAACTG 






FK832766 






TGAGGGAACTCCCTGTTCAC 


Down 
syndrome 


DSCAM 


Coptotermes 
formosanus 


Downsyndrome cell 
adhesion molecule 


immunglobulin- 
like fibronectin. 


CCGTACCATTGACGGAGAGT 


cell adhesion 
molecule 




FK835436 


Harpegnathos saltator 
(EFN81 146.1} 


cell recognition 
and adhesion. 


CACTTGTGTTGCCCGTAATG 


axon guidance 
receptor in insects 


Ficolin-2 
precursor 


FICO-2 


Coptotermes 
formosanus 


Flcolin-2 precursor, Pediculus 
humanus corporis, 
(XP_002430827.1) 


Antigen 
recognition, 
carbohydrate binding, 


CTCCGTCGGACTCTCATCTC 






FK836944 


scabrous protein-like, Tribollum 
castaneum, (XP_972571.1} 


Enhances 
phagocytosis and 
opsonization 


GCGCTATAACCTGGGTCTCA 


Lectin-related 
hemo lymph 


LECT-like 


Coptotermes 
formosanus 


Hemolymph lipopolysaccharide- 
binding protein, Periplaneta 
americana (P26305.1) 


£ coli 

Lipopolysaccharide 
binding. 


CTTCCCAGTGTCGAAAGCTC 


Lipospolysacch 
aride binding 
protein 




FK834461 


Lectin-related regenectin, 
P. americana (BAA82265.1) 


Tissue 

regeneration 


TGGGTACCCCATAGGAGATG 


C-type lectin 


CLECT 


Coptotermes 
formosanus 


conserved hypothetical 
C-type lectin, P. tiumanus 
corporis (XP_002424769.1) 


carbohydrate- 
binding, pathogen 
recognition 


AGGGTGAAAGAGGAGCAACA 






FK835521 




innate immunity 


TAAAAATCCCAACAACCCCA 


Gram-negative 
binding protein 2 


GNBP-2 


Drepanotermes 
rubriceps 


gram negative bacteria binding 
protein 2 from different termite 
species 


insect pathogen 
recognition protein 


CTTCCCAGTGTCGAAAGCTC 






DQ058934 


(e.g., AAZ0893.1-AAZ08504.1, 
AEK64801.1, ADJ19004.1) 


beta(1,3)- 
glucanase effector 


TGGGTACCCCATAGGAGATG 


Reference 
genes 


Cytoplasmic 
heat shock 


HSP 


Coptotermes 
formosanus 




ATP and 

nucleotide binding 


AGCCTTGGCCACAACAGTGCAA 



protein 70 
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Table 1. Cont. 





Gene 


Abbrev. 


Species 
of Origin 


Top Matches in GenBank, 
Species of Origin 


Putative 
Function 


Primer 

Sequence (5'-3') 






(Accession No.) 


(Accession No. for Nucleotide 
or Peptide Sequence) 










FK835495 






GGAGCAGGAGCAGGACCGACT 


NADH 

Dehydrogenase 
subunit 4 


NADH 


Coptotermes 
formosanus 




Mitochondrial 
electron, sodium 
ion and 


ACGAAGCAACCCATAACCACCAAGC 






FK833785 




proton transport, 
ubiquinone activity 


GGGCTCATGTTGAGGCTCCTGTT 


Elongation 
factor! -alpha 


EFA 


Coptotermes 
formosanus 




GTP binding. 


GGCGGTCACATTTCTCCTTA 






FK834645 




regulation of 

translational 

elongation 


CTGAACCACCCTGGTCAGAT 



doi:1 0.1 371/journal.pone.0105582.t001 



0.0001) in E. coli/P. lermilis injected termite.s than in 0.5% PB fed 
termites. MatE was only induced (F= 32.37, P<0.001) in the E. 
coli/P. termitis treatment group, while gene expression after PB 
treatment was not different from that of the controls. 

In contrast, GST was exclusively induced after 0.5% PB 
treatment (F = 3.53, P = 0.03 12), while gene expression did not 
differ from controls after E. coli/P. termitis treatment. Although 
induced by both treatments, relative expression of LECT-like, 
CLECT and FICO-2 was significandy higher (F = 33.03, P< 
0.001; F = 61.72, P<0.0001; and F=23.78, P<0.001, respective- 
ly) for termites fed with 0.5% PB than for termites injected with E. 
coli/P. tennitis. 

4. Effects of caste 

Tests for differential expression showed that caste (worker or 
soldier) had a significant (P<0.05) effect on levels of expression in 
five of the eight inducible target genes (CYP15A1, MatE, GST, 
GNBP, and CLEG, Figure 1, Table SI). The gene expressing 
CYP15A1 was induced in both workers and soldiers, with a 
significant interaction (P = 0.032) of treatment and caste, due to 
the significantly higher (F= 34.79, P<0.0001) expression in 
workers as compared to in soldiers. Gonstitutive expression of 
CYP15A1 was also significantly increased (F= 213.65, P<0.001, 
PROG GLM) for untreated workers as compared to untreated 
soldiers. 

Gonstitutive expression of MatE was significandy higher 
(F = 40.16, P<0.0001) in untreated workers than in untreated 
soldiers as well. Furthermore, MatE was only induced in workers 
by septic injury. The Gram-negative bacterial binding protein 
(GNBP-2) was, similarly, almost exclusively induced in workers in 
both treatments (F= 121.23, P<0.001). GST was also exclusively 
induced in workers, but only after challenge with 0.5"/o PB. The C- 
type lectin gene (CLEGT), was the only target gene which showed 
significandy higher (F = 47.53, P<0.001) expression in soldiers 
than in workers. 

Gaste had a marginal effect on the inducibUity of Ficolin-2 
(FICO-2) and chitin binding protein (GBP). These genes were 
induced both in workers and in soldiers, but FIGO-2 was 
expressed marginally higher in the worker caste (F = 7.31, 
P = 0.058), while GBP's expression was marginally biased towards 
soldiers (F = 8.10, P = 0.049). The second lectin related gene 
(LECT-like) also was induced in both castes, but showed no 



significant difference in expression levels as an effect of caste 
(F = 0.16, P = 0.692). 

Discussion 

1. Effect of xenobiotic exposure and septic injury on 
expression of target genes 

Constitutive expression of all the investigated genes in EST 
workers and soldiers was low; however, eight out of the ten genes 
were inducible by one or both treatment types in one or both 
castes. The initial categorization based on the putative function of 
the ten target genes into genes associated with detoxification/ 
xenobiotic metabolism and septic injury/ immunity response was 
somewhat supported by the induction experiments using pheno- 
barbital (PB) treatments and septic injury with bacteria, however, 
there were also some unexpected results (Figure 2). 

The expression of the glutathione S-transferase (GST) and one 
of the cytochrome P450 genes (GYP15A1) as expected was 
upregulated in response to exposure with the xenobiotic PB. The 
barbiturate PB is used to treat epilepsy in humans and acts on 
GABA-gated chloride channels in the central nervous system to 
reduce excitatory synaptic response similar in activity to some 
termiticides, i.e., FipronU [33]. 

Response elements and nuclear receptors responsible for 
induction of many cytochrome P450s and GSTs have been 
described in insects and shown to respond to several groups of 
xenobiotics, including PB [34]. In insects the GSTs are known to 
metabolize plant allelochemicals and insecticides, and high levels 
of GSTs are related to insecticide resistance. Exposure to PB has 
been previously shown to induce GST activity in cockroaches [35], 
Drosophila [36,37,38], ^efi(« aegypti [39], among others (reviewed 
in [40]). Upregulation of the expression of the GST gene in FST 
workers upon exposure to PB was therefore expected. The lack of 
upregulation in soldiers and the effect of caste on differential gene 
expression are discussed below. 

Although, cytochrome P450 enzymes are largely ubiquitous and 
have common catalytic mechanisms, they are translated from a 
large family of genes with high sequence diversity, which explains 
their numerous functions. In insects, cytochrome P450s are known 
to play a role in detoxification of insecticides and xenobiotics 
(reviewed in [40]), but also in regulation of caste development via 
metabolism of hormones, like the juvenile hormoneJH [41,42,43]. 
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A* .A* ^^5.-* ^"t* 



Figure 1. Constitutive and induced expression of target gene amplicons. The y-axes denote relative expression in R units averaged over 
three colonies (note that the ranges of y-values vary across graphs). The x-axes denote caste (W = worker, S = soldier) and treatment type (+PB = fed 
with 0.5% Phenobarbital, +B = septic injury with £ coli/P. termitis and -PB/-B = untreated control showing constitutive expression). 'IVlissing' columns 
have a relative expression of<0.01. + marks significant induction. There was no induced expression observed for either CYP450 or DSCAM (not 
shown). 

doi:10.1371/journal.pone.0105582.g001 



The closest matches for the gene product of one of the CYP450s 
in our study were a CYP15A1 of the cockroach Diploptera 
punctata (AAS13464.1, [41]) and a CYP15A1 pseudogene found 
in R. flavipes (ACN93795.1, [43]) with an identity of ,51-54%, 



which lies within the range of cross-species alignments of 
translated CYP15s (41,5-72.5%, [43]). The product of CYP15A1 
epoxidizes methyl farnesoate to juvenile hormone III in cock- 
roaches [41,44] and a similar function in JH biosynthesis was 
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suggested in termites [43,45]. Tarver et al. [43] sliowed tliat 
CYP15F1, wliicli is closely related to the truncated CYP15A1 
pseudogene, is involved in JH dependent caste differentiation. 
However, this study [43] did not investigate whether these two 
genes were ffinctionaUy related. 

Unlike other prominent cytochrome P450s, which metabolize a 
variety of pesticides due to their broad substrate and product 
specificity, CYP15A1 was shown to be highly specific in the 
silkworm [41]. The CYP15A1 gene in our study was the only 
target gene that was constitutively expressed and induced by both 
treatment types in workers and soldiers. The fijndamental role in 
JH biosynthesis and caste regulation of similar genes in termites 
[43] may explain the level of constitutive expression observed in 
our study. 

The inducibility by PB and septic injury, however, suggests that 
the CYP15A1 gene product has expanded functionality in terms of 
xenobiotic metabolism and immunity beyond its primary function 
as epoxidase catalyzing the last step in JH biosynthesis. While PB 
inducibility of CYP450s of the families 4, 6, 9, 12 and 28 in insects 
is well known, the only inducers for C YP 1 5 family members so far 
reported were JH and terpene components from soldier head 
extracts suspected to act as caste regulatory primer pheromones 
[43]. Our study is the first to demonstrate inducibility of a member 
of the CYP15 family by a xenobiotic. However, the involvement of 



the CYP15 enzyme in actual xenobiotic metabolism and 
detoxification has to be confirmed in enzyme and bioassays. 
Similar to CYP15A1, another enzyme involved inJH metabolism 
{JH III epoxide hydrolase responsible for JH degradation) has been 
shown to be upregulated after PB exposure in several insect species 
([46], and references therein). 

A link between JH levels and innate immunity has been 
discovered in the silkworm Bonibus mori [47], when JH treatment 
acted as an immune activator and increased the transcription 
products of six AMPs. We suggest that septic injury in our study 
upregulated CYP15A1 expression in termite workers and soldiers, 
which results in JH biosynthesis to provide increased JH levels as 
an immune activator. This route would constitute an immune 
pathway different but complementary to the ToU- and Imd 
pathways. 

Overexpression of CYP450s and GSTs is often associated with 
insecticide resistance ([37], and references therein). However, the 
investigated CYP450s and GST genes were not constitutively 
expressed at elevated levels in Formosan termites and induction 
with PB led to a rather moderate spike in gene expression (less 
than 4 fold), which is in the range of susceptible insects [48]. 
Insecticide resistant strains usually show a much higher fold- 
change [49]. Although there is variability in tolerance to 





[ W 




[ Xenobiotic 


+ - ^ 
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Xenobiotic 
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Figure 2. Graphic summary of caste- and treatment dependent inducibility and function of eigKit detoxification and immune- 
related genes of the Formosan subterranean termite. Treatment consisted of exposure to a xenobiotic substance (0.5% phenobarbital) or 
septic injury via bacteria injection. Relative levels of caste (workers vs. soldiers) dependent induced gene expression are indicated with +++ in red (>4 
fold expression change compared to control), ++ in orange (2-4 fold expression change), + in yellow (<2 fold expression change) or - in grey (no 
significant induction). Genes are arranged from left to right according their putative function ranging from detoxification to immunity and from top 
to bottom according to their caste dependent expression. 
doi:1 0.1 371/journal.pone.01 05582.g002 
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insecticides among termite colonies and species, to date, no 
resistance or treatment failure has been recorded. 

The second cytochrome P450 gene in the study (CYP450) did 
not respond to treatment at all. CYP450 in our study was most 
closely related (48 'Xo identity) to shade (CYP314al) from the 
Halloween gene family, an ecdysone 2()-monooxygenase ecdysone 
hydroxylase involved in 2()-hydr{)ecdysone (molting hormone) 
synthesis in D. melanogasler [50]. GYP450 also showed consid- 
erable similarity to CYP12A2 from the house fly Musca domestica. 
A similar gene (CYP12A1) in the house fly was shown to be 
inducible by PB and to mc-tab()liz(- different insecticides [51]. 
Thus, we initially hypothesized that this gene would be most likely 
up-regulated in response to PB treatment and possibly involved in 
xenobiotic metabolism. However, this gene apparently lacks the 
response elements and/or transcription factors needed for PB 
inducibUity [34]. Genome-wide studies in other insects showed 
that only a third of GYP genes are inducible by xenobiotics [52]. 

Similar to GST and GYP450 genes, we initially expected the 
gene for MatE to be inducible by the xenobiotic PB as MatE 
proteins have been described as multi drug resistance efflux 
proteins. They function as transmembrane pumps and are 
involved in the excretion of xenobiotics in vertebrates, detoxifica- 
tion of metals and secondary metabohtes in plants [53] and 
mediating resistance to antibiotics in bacteria [54,55]. It was 
somewhat unexpected that MatE expression was upregulated by 
septic injur\' and not by xenobiotic challenge in our study. The 
only evidence to date that MatE might be involved in antimicro- 
bial defense is found in plant studies. Nawrath et al. [56] 
investigated a MatE gene homologue in Arahidopsis plants: 
expression was strongly induced by inoculation with an avirulent 
strain of Pseudomonas syringae and the authors concluded that 
this gene was likely involved in antimicrobial plant defense via the 
salicylic acid-dependent signaling pathway. 

To our knowledge, MatE proteins and their function in insects 
have not yet been investigated. It is possible that the MatE gene in 
our study was directly involved in immune response in termites as 
indicated by the induced expression in EST workers injected with 
E. coli and P. termitis. However, it is m()r(; likely that MatE was 
upregulated as a consequence of the production of cytotoxic 
metabolites resulting from the activation of immune pathways by 
bacteria injection. For example, septic injur)' in insects activates 
the PPO pathway, which produces melanin, a phenohc biopoly- 
mer. Some melanin precursors are cytotoxic, which is advanta- 
geous for killing opportunistic invaders at the wound site [57]. 
However, metabolites might be inadvertentiy cytotoxic to the 
insect itself, especially since the PPO cascade is taking place in the 
open circulatory system and thus byproducts may accumulate in 
the hemolymph. MatE proteins play an important role in the 
removal of metabolic waste products and are believed to be a 
transport system for phenolic compounds [58]. Excessive melanin 
or its intermediaries from the PPO pathway might thus be 
extruded or sequestered by MatE in termites. 

We initially assumed that the cliitin binding protein (GBP) 
belongs into the category of immune related genes. The GBP in 
our study matched closest (31%) to a putative peritrophic 
membrane protein in a scarab beetle {Holotrichia oblita.), which 
contains carbohydrate and chitin binding domains. Invertebrates, 
such as shrimp, possess AMPs with chitin-binding capacity 
thought to be involved in bacteria phagocytosis, and perhaps 
antifungal activity and wound healing via chitin assembly [59]. 
Upregulation of GBP upon septic injury in workers and soldiers of 
Formosan termites suggests a similar involvement in immunity 
and/ or wound healing via restoration of the cuticle. 



However, PB treatment also induced GBP expression albeit at 
lower expression levels than bacteria challenge. A genome-wide 
analysis in the flour beetle, Triholium castaneum, revealed that 
some proteins with peritrophin A-t)'pe chitin-binding domains 
were expressed in cuticle-forming tissues while others were 
associated with the peritrophic membrane and expressed in the 
midgut during feeding [60]. The peritrophic membrane is a 
barrier against pathogens and a sink for toxic substances [61], 
which would explain upregulation of certain peritrophins after 
ingestion of a xenol)iotic like PB. 

Similarly, the Ficolin 2 precursor (FIGO-2) and the two lectins 
(GLEGT, Lectin-Uke) in our study were assumed to be involved in 
immune defense enhancing phagocytosis and opsonization [62]. 
Ficolins are fibrinogen related proteins that act as PRRs by virtue 
of binding to carbohydrates on the surface of gram-negative and 
gram-positive bacteria [63]. This leads -at least in vertebrates- to 
the activation of the lectin complement pathway resulting in 
bacteria lysis. In insects, ficolins are not that well studied but the 
fact that their fibrinogen-like domain is conserved throughout the 
animal kingdom suggests that they are involved in ancestral forms 
of the innate immune system in invertebrates [64-66]. Genome- 
wide analyses and microarrays in Drosophila spp [15,66] 
confirmed that fibrinogen-related domains are involved in 
immune response in insects. The Ficolin precursor in our study 
showed high similarity (78%) to a putative Ficolin-2 precursor of 
the human louse {Pediculus humanus corporis) and scabrous 
protein of T. castaneum (72%) and was upregulated by septic 
injury. 

Both lectin related gene products in our study belong to the G- 
type lectin superfamily of calcium dependent carbohydrate- 
recognition proteins. G-type lectins play a variety of roles in 
adaptive and innate immunity as PRR, cell adhesion molecules 
and activators of the complement pathway [67]. Lectins in the 
hemolymph of cockroaches and lepidopterans participate in innate 
immune response via phagocytosis, encapsulation, nodulation and 
activation of the PPO pathway ([9], and references therein). 
Immunelectins vary in their specificity towards microbial epitopes 
as some show induction after Gram-negative bacteria injection 
only, while others bind to Gram-negative and Gram-positive 
bacteria as well as fungi [9] . 

The product of the LEGT-like gene in our study showed 32% 
identity to both a hemolymph lipopolysaccharide (LPS)-binding 
protein [68] and regenectin [69] from the cockroach Periplaneta 
americana. The LPS binding protein had affinity for lipopolysac- 
charides on the cell membrane of some Gram-negative bacteria, 
but not others, suggesting high selectivity to certain oligosaccha- 
ride structures and thus specific immune recognition of selective 
pathogens [68]. Regenectin is involved in muscle formation and 
leg regeneration in cockroaches, but in contrast to LPS-binding 
protein, the transcription of its gene is not enhanced by bacterial 
infection [69]. Transcription of the LEGT-like gene was 
upregulated in EST soldiers upon septic injury, which most likely 
reflects bacteria binding function and/ or participation in tissue 
regeneration for wound heafing. 

The product of the second lectin related gene (GLEGT) in our 
study was highly related to a conserved hypothetical protein in the 
human body louse (90% identity) and G-type lectin-like precursors 
from several insect species such as the flour beetie several 
mosquito species (S86%), Drosophila spp. (a86%), and social 
insects, i. e., ants, honey and bumble bees (S88%). GLEGT 
belongs thus to the immunelectins acting as PRR [9] that activate 
the PPO pathway. C-type lectins are also involved in sclerotiza- 
tion, i.e., exoskeleton hardening in arthropods [70]. Sclerotization 
is catalyzed by phenoloxidases and thus linked to activation of the 
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PPO pathway. Furthermore, the calcium binding properties of C- 
type lectins are assumed to contribute to biomineralization of the 
cuticle by depositing calcium salts [70]. Thus, it was not surprising 
that CLECT was upregulated after septic injury similar to the 
LECT-like gene. 

We were surprised, however, that the expression of the Ficolin 
precursor and both lectins was also induced by PB treatment. The 
majority of antiepUeptic drugs such as PB show immunosuppres- 
sive effects in humans; however, under certain conditions, they can 
also stimulate the immune .system [71]. Microarray analysis of 
gene expression revealed that Drosophila flies execute massive 
transcriptional response to PB treatment [72]. Not only were the 
detoxificication genes upregulated, but also genes involved in 
several metabolic and stress response pathways were affected. A 
nuclear receptor (DHR96) played a central role in regulating the 
response to PB. Middha and Wang [66] subsequently .showed that 
some fibrinogen-related proteins are regulated by DHR96, thus 
Unking expression of fibrinogen-like genes to PB exposure. This 
suggests a broader than originally assumed network of functions 
beyond innate immunity mediated by fibrinogen-like proteins in 
insects. 

Gram-negative binding proteins are known in insects to sense 
pathogen invasion and trigger immune responses, i.e., the Toll and 
Imd pathways [73], and references therein) as well as the PPO 
cascade [74] . In termites, GNBP-2 has a double function as PRR 
that recognizes LPS (like lectins) and beta-l,3-glucans and as 
antifungal effector due to its beta-l,3-glucanase activity [21]. 
GNBP-2 is expressed in the salivary gland and by hemocytes and 
secreted via allogroming onto the cuticle and incorporated into 
nest and tunnel material presumably to suppress growth of fungal 
pathogens [21,73]. 

Prior to our study, the role of GNBP in termites has been 
explored mostly in terms of antifungal but not antibacterial activity 
[73,75,76]. Our study showed involvement of GNBP-2 in defense 
against bacterial septic injury deduced from the fact that 
expression in workers doubled after bacteria injection. This result 
supports Bulmer et al.'s [21] in vitro binding assays that showed 
that GNBP-2 had affinity to Gram-negative bacteria. 

While the upregulation of GNBP-2 upon immune challenge was 
expected, it was surprising that feeding on PB also induced C!NBP- 
2 expression in workers. There is no straightforward explanation, 
but we may speculate that workers increase grooming of colony 
mates showing signs of PB exposure, which might indirectly lead to 
increased production of GNBP-2 in the salivary gland as 
protective secretion [21,73]. Furthermore, PB acts as universal 
inducer for many genes in insects, not only CYP450s and GSTs 
but also for juvenile hormone epoxide hydrolases and for genes 
involved in carbohydrate metabolism, such as amylases [37]. 
GNBP-2 has a glucanase site that performs hydrolysis on 
glucosidic bonds and thus is a hydrolase involved in carbohydrate 
metabolisms. Like the increased transcription of other PRR 
immune genes upon xenobiotic exposure (see lectins and ficolin 
above) the regulatory and functional link between detoxification 
and immunity needs to be further investigated. 

We initially hypothesized that the Down syndrome cell adhesion 
molecule (Dscam) was immune related. The gene in our study 
showed several putative conserved domains of the Immuneglobu- 
lin superfamUy and fibronectin Type III similar to those in 
Hymenoptera and Diptera [77], which play a role in cell adhesion 
and wound healing in vertebrates. Insect studies showed that the 
Dscam receptor is involved in neural wiring and innate immunity 
[78]. Reduction of Dscam expression using mutants or RNA 
interference impaired phagocytosis of bacteria presumably due to 
reduced binding and opsonization of bacteria [79,80] . 



Surprisingly, Dscam was not induced by septic injury (or by PB 
treatment) in our study. Dscam in the insect genome has many 
alternative exons and encodes tens of thousands of isoforms 
resulting from alternative splicing [77,78]. Watson et al. [80] 
showed that some Dscam isoforms bound to E. c.oli, while others 
did not. There might be bacteria binding and non-binding 
isoforms in general, and/or binding is somewhat specific to 
epitopes of non-resident bacteria species, e. g., pathogens. The 
bacteria in our assay were native symbionts cultured from the 
termite gut, so the Dscam in our study might not be triggered. 

2. Gene expression levels are mostly biased in favor of 
workers 

Expression of five of the eight inducible target genes (GYP 1 5 A 1 , 
GST, MatE, GNBP-2, GBP and CLECT) was significandy and 
that of two additional genes (FICO-2, CBP) was marginally 
influenced by caste. Constitutive and/ or induced gene expression 
was higher in workers in most of those differentially expressed 
genes except for CLECT, which was predominantly expressed in 
the soldier caste. In fact, three of the target genes (GST, MatE, 
and GNBP-2) were almost exclusively induced in workers. 

In our study, higher expression of the "classic" detoxification 
genes such as members of the CYP450s (CYP15A1) and GST in 
the worker caste supports previous findings of higher enzyme levels 
in workers and higher tolerance against insecticides. Gatti and 
Henderson [6] and Gatti et al. [81] showed higher tolerance 
against organophosphates in workers attributed to a higher 
acetylcholinesterase activity in workers. Valles and Woodson [8] 
discuss their unpublished data which showed that soldiers of C. 
formosanus have only a fraction of detoxification enzyme activity 
compared to workers. Nevertheless, in a study of Osbrink et al. [7] 
soldiers of C. formosanus and R. virginicus in several colonies 
were more tolerant than workers to the majority of insecticides 
tested. This obvious variance among studies, species, enzymes and 
insecticides underlines the need for continued investigation of the 
relationship among gene expression, inducibUity, metabolic 
enzyme production and activity, and insecticide susceptibility. 

The major immune-related genes (GNBP-2, FICO-2) were 
induced and upregulated at higher levels after septic injury in 
workers than in soldiers. Soldiers might not rely primarily on the 
production of their own immune peptides for protection against 
pathogens. Soldiers have defense secretions, which in some species 
may have antifungal properties [23] and a higher degree of cuticle 
sclerotization than workers, which may aid as a mechanical 
barrier. Soldiers also ma)' lia\c antibacterial and/or antifungal 
activity via chitin assembly, wound healing and systemic protec- 
tion related to increased CBP expression, which has been 
suggested for other highly sclerotized arthropods [59]. Sclerotiza- 
tion results from melanization and melanin plays an important 
role in innate immunity via wound healing and encapsulation of 
invaders. Melanin production during the PPO pathway also leads 
to toxic metabolites, such as quinones that could kill invaders [82] . 
Thus soldiers might have more non-protein related mechanical 
and biochemical defenses in form of sclerotization, melanin and its 
byproducts in addition to antimicrobial gland secretions. Soldiers 
could also receive passive immunization via secretions from 
workers. It has been shown that GNBP-2 is produced in the 
workers' salivary gland and secreted to coat nest material and the 
cuticle [21] and that social transfer of immunity can occur [83]. 
Workers might overproduce GNBP-2 and other immune peptides 
and transfer them to soldiers via grooming. 

CLECT was the only gene that was expressed at significandy 
higher levels in soldiers. C-type lectins are involved in mineral- 
ization and sclerotization of the cuticle [70] . Soldiers are more 
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sclerotized than workers and thus upregulation of CLECT 
expression might be naturally facilitated and high levels of 
CLECT can be produced on demand. Although soldiers are 
likely to frequently get injured during colony defense, they did not 
upregulate the other immune proteins investigated in this study 
(GNBP-2 and Ficolin) to the same degree as workers. Conse- 
quently, the soldier caste might primarily rely on lectins, which 
they use naturally for exoskeleton hardening to support their 
innate immunity and complement passive transfer of immune 
proteins from workers secretions. Research into caste specific 
differences, concerning not only the level of immune defense but 
also the mechanisms and pathways used, might be of interest in 
future studies on social insects. The present study provides the 
foundation for this endeavor. 

Materials and Methods 

1. Sample collection 

Worker and soldier FSTs were collected in May 2009 from 
three field colonies in New Orleans, Louisiana, USA, an urban 
habitat in which this invasive species thrives: Colony 1 - City Park, 
adjacent to Tad Gromley Stadium at 1 Palm Drive; colony 2 - 
Louis Armstrong Park at 90 1 North Rampart Street and colony 3 - 
French Quarter at 300 Canal Street. No specific permissions were 
required for accessing these locations for sampling activities and no 
endangered or protected species were involved in the studies. 
Voucher specimens are deposited in the Louisiana State Arthro- 
pod Museum at the LSU Agricultural Center. Based on previous 
studies of FST colony size and flight distance [84—86] , samples 
from these locations represented three independent colonies from 
the New Orleans population. Termites (~500 workers and 500 
soldiers) were collected from untreated inground monitoring 
station (Sentricon, Dow AgroSciences LLC: Indianapolis, IN) at 
each location and then transferred to plastic tubs (50x30x20 cm) 
lined with moist corrugated cardboard and maintained overnight 
in the laboratory at 25 °C and 30% relative humidity (InteUus 
Environmental Controller, Percival: Perry, lA). 

2. Determining sublethal treatment conditions for septic 
injury and phenobarbital challenge 

Since bacteria invasion is generally associated with injury, and 
the responses to injury and infection are likely to be intricately 
connected, we decided to test the response to septic injury. To 
determine the appropriate sublethal treatment for inducing gene 
expression related to septic injury non-pathogenic bacteria, that 
were previously cultured from the guts of Formosan termite 
workers (gram-negative Escherichia coli [87] and gram-positive 
Pilibacter termitis [88] ATCC BAA- 1030), were revived from 
glycerol stocks. E. coli was cultured under aerobic conditions in 
1 mL of BHI broth for 1 day on a shaking incubator at 30°C, 
while P. termitis was grown anaerobically, using H2/CO2 BBL 
GasPaks (Becton Dickinson and Company: Franklin Lakes, NJ), in 
1 mL of BHI broth at 30°C for 3 days. Sterile 28 gauge needles 
(Comfort Point: Nanaimo, British Columbia) were each soaked for 
approximately 60 seconds in a mixture of equal volumes of 
exponential growth phase E. coli and P. termitis (200 |a,L 
supernatant of each, placed in a 1.5 mL micro-centrifuge tube). 
Twenty workers and 20 soldiers from colony 1 were cold 
immobilized on a thermoelectric laboratory chill table (Bioquip: 
Gardenia, CA) and placed with their ventral surface exposed on a 
dissecting microscope stage (Leica MZ7.5 Stereomicroscope, 
Meyer Instruments: Houston, Texas). The abdomen was swabbed 
with 70% ethanol before being pierced through the intersegmental 
membrane of the 4* and 5'^ segments with the septic needle [18]. 



Following inoculation with bacteria, individual termites were 
isolated in 47 mm Petri dishes and kept at 25°C and 30% relative 
humidity (InteUus Environmental Controller, Percival: Perry, lA) 
for 48 hours. Twenty workers and 20 soldiers, which were not 
inoculated with bacteria, were held individually in separate dishes 
under the same conditions as controls. Mortality was determined 
as the inability for the termite to right itself after 30 seconds and 
was assessed for treated and untreated termites at 12, 24 and 
48 hrs. Since no mortality was recorded at 24 hrs but 50% and 
70% mortality in workers and soldiers, respectively, was recorded 
at 48 hrs, the 24 hr treatment was chosen as the sublethal 
treatment condition for target gene induction (Table 2). The 24 hr 
sublethal exposure time was consistent with rapid induction of 
innate immunity in insects (i.e., peptide synthesis in Drosophila 
[89] and in P. spiniger [18]). 

For the xenobiotic challenge we chose phenobarbital (PB, 
Sigma-Aldrich: St. Louis, MO), a universal inducer of a broad 
variety of detoxification genes in insects [36,40]. Most studies test 
reactions to xenobiotics by topical application or ingestion. Since 
those would be the natural pathways of exposure in termites we 
used droplet feeding, which likely leads to ingestion and topical 
exposure via contact with the droplets and possibly transfer by 
grooming. To determine sublethal dose and time for PB 
treatments, 20 workers and 20 soldiers from colony 1 were each 
placed individually in a 47 mm Petri dish (MUlipore: BiUerica, 
MA). Five 2 |J.L droplets of 1% PB were pipetted directly onto the 
inside surface of each Petri dish. Droplets of 1% PB also were 
added to 20 Petri dishes which did not contain a termite as an 
evaporation control. All 60 dishes were sealed with parafilm and 
kept at 25°C and 30% relative humidity (InteUus Environmental 
Controller, Percival: Perr^-, lA) for 48 hrs. The PB-infused water 
droplets were monitored at 12 hr interx^als. Level of evaporation 
was negligible over the 48 hr incubation period, and the PB 
droplets were consumed by workers and soldiers at approximately 
equal rate, confirming uptake. 

We then used the droplet feeding assay to determine the 
maximum sublethal time and dose of PB treatment. Phenobarbital 
was serially dUuted with distiUed water (1%, 0.5%, 0.25% and 
0.125%), and each concentration of PB was fed to 20 workers and 
20 soldiers in 2 nl droplets as described above. Controls were fed 
distilled water only. Mortality was assessed as described above. 
Based on the results of this experiment (Table 2), the 0.5% PB 
treatment at 24 hrs was chosen as the sublethal treatment 
condition for target gene induction in both workers and in 
soldiers, which lies within the published time frame recommended 
for measuring PB- induced gene expression in insects ranging from 
4 hrs in Drosophila [38] to 5 days in cockroaches [35]. 

3. Experimental design for septic injury and PB challenge 

FoUowing the established sublethal treatment conditions, 20 
workers and 20 soldiers from each of the three FST colonies were 
fed 0.5% PB, 20 workers and 20 soldiers per colony were 
inoculated with the E. coli/P. termitis mixture and 40 workers as 
weU as 40 soldiers per colony were left untreated (controls). 
Termites were kept in individual Petri dishes at conditions 
described above and flash-frozen at — 80°C (Ultra low VIP, Sanyo 
Scientific: Wood Dale, IL) after 24 hrs. From each of the three 
colonies (biological replicates), five individuals of each caste 
(worker/ soldier) and treatment (0.5% PB, septic injury with E. 
coli/P. termitis untreated control) were randomly chosen as 
technical replicates for quantitative real time PGR (qRT-PCR) 
reactions. 
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Table 2. Mortality of FST workers and soldiers to determine sublethal treatment conditions. 



Time (h) 


Treatment 


Mortality (%) 








Workers 


Soldiers 


48 


1% PB 


100 


100 


24 


1% PB 


70 


70 


12 


1% PB 


50 


70 


48 


0.5% PB 


50 


50 


24 


0.5% PB 


0 


0 


12 


0.5% PB 


0 


0 


48 


0.25% PB 


0 


0 


24 


0.25% PB 


0 


0 


12 


0.25% PB 


0 


0 


48 


0.125% PB 


0 


0 


24 


0.125% PB 


0 


0 


12 


0.125% PB 


0 


0 


48 


E.coli/P. termitis 


50 


70 


24 


E.coli/P. termitis 


0 


0 


12 


E.coli/P. termitis 


0 


0 


48 


Control 


0 


0 


24 


Control 


0 


0 


12 


Control 


0 


0 



Termites were fed with specified concentrations of serially diluted phenobarbital (PB), injected with a 50:50 mixture of f. coli and P. termitis or left untreated (Control). 
Sublethal treatment conditions {maximum dose and/or exposure time prior to observed mortality) for target gene induction are in bold print. 
doi:1 0.1 371 /journal.pone.Ol 05582.t002 



4. RNA isolation and cDNA synthesis 

Total RNA was extracted from whole bodies of individual 
termites using the RNA Purification Microkit (QIAGEN, 
Valencia, CA). RNA concentrations were assessed at the 260/ 
280 nm ratio of absorbance using a ND-1000 spectrophotometer 
(Nanodrop Technologies: Wilmington, DE). RNA samples were 
normalized to 10 ng per |a.l by adding RNAse free water and were 
then used as template for complimentary DNA (cDNA) synthesis 
via the iScript Kit (Biorad: Hercules, CA, USA) using a PTC- 
200 DNA Engine gradient thermocycler (MJ Research: San 
Francisco, California). The cDNA was normalized to 1,000 ng 
per |j,l and used as template for qRT-PCR. 

5. Primer design for target and reference genes 

We selected ten target genes putatively involved in immunity 
and xenobiotic metabolism from our Expressed Sequence Tag 
(EST) library of the FST [29], the literature and the NCBI 
database. Two ESTs matching to genes from Cytochrome P450 
families [CYP15A1 (FK835449), CYP450 (FK833823)], and ESTs 
encoding for glutathione S-transferase [GST (DC239424)] and the 
multi antimicrobial extrusion protein [MatE (FK833694)] were 
assumed to be mainly involved in xenobiotic metabolism. Two 
ESTs encoding for lectin-related peptides [CLECT (FK835521), 
LECT-like (FK834461)], and those encoding for a Ficolin-2 
precursor [FICO-2, (FK836944)], the Down syndrome ceU 
adhesion molecule [Dscam (FK835436)], a chitin binding protein 
[CBP (FK832766)] and the gram-negative binding protein 2 
[GNBP-2 (DQP58934)] are likely to play a role in pathogen 
recognition due to their putative cell adhesion function and were 
thus assumed to be immune related. 

Primers for GST and GNBP-2 were already developed for a 
previous study that examined differential gene expression among 



FST queens and precopulatory females [29]. Three primer pairs 
were designed for each of the remaining seven target genes 
(PRIMER-BLAST software suite, NCBI) to have between 17 and 
28 bases, with -50% to 60% GC content. 

The EST library for FST [29] was first screened for reference 
genes used in previous studies (e.g., [90,91]) for quantifying gene 
expression, i.e., Glyceraldehyde-3-phosphate dehydrogenase 
(FK834351), TATA-binding protein (FK835420), Elongation 
factor- 1 -alpha (FK834645) and Ubiquitin specific peptidase 
(FK831079). Next, the library was screened for genes most likely 
to be involved in cell maintenance and/ or basal cell metabolism, 
i.e., "housekeeping" genes. Cytoplasmic heat shock protein 
(FK835495) and NADH dehydrogenase (FK833785) were select- 
ed, bringing the total number of putative reference genes to six. 
Table S2 lists the two to three primer pairs designed for each of 
the reference genes, along with their melting points, sequences and 
GC content. 

Each primer pair designed for target genes was tested using 
qRT-PCR with cDNA from three untreated, randomly chosen 
workers from colony 1 as template. Primer pairs designed for 
reference genes were tested with cDNA template synthesized from 
three colony 1 workers and soldiers that were randomly chosen 
from each of the 0.5% PB treatment group, the E. coli/P. termitis 
treatment group and the untreated control group. 

Amplifications of target and reference genes were performed 
using the IQ_5 iCycler (Biorad: Hercules, CA, USA) in reactions 
each containing 5 )J.l cDNA, 25 |a.l SYBR-green Supermrx 
(BioRad), 150 nm forward and reverse primer for each candidate 
or reference gene (see below), and 1 7 [il H^O for a total volume of 
50 |xl. The thermal protocol consisted of one cycle of 95°C for 
3 mill, followed by 40 cycles of 1) 95°C for 10 sec and 2) 58°C for 
30 sec. Relative fluorescence units and Ct values (fractional cycle 
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at which amplification reaches a detection threshold) were 
recorded at the end of each cycle by the iCycler software. 

For each target gene, 1 of the 3 primer pairs was selected based 
on the following criteria [92-94]: 1) consistency across technical 
replicates (different individuals run on different plates), 2) target 
amplicon at 50 to 100 bp, 3) melting curve with a single peak and 
4) reaction efficiency (E) Sl.8. Primers for reference genes needed 
to pass all of the same criteria while also being consistent across 
treatment types (0.5% PB and E. coli/P. termitis, control) and 
among castes (workers and soldiers). Table S2 lists the ten target 
and three reference genes with their corresponding primer pairs 
which passed all of the above criteria. 

6. Relative expression and tests for differential expression 

LinRegPCR software [93] was used to calculate the individual 

efficiencies (E) for reactions, which were then averaged across the 
five technical replicates for each gene and treatment type. Cx 
values given by the iCycler qRT-PCR software were similarly 
averaged across technical replicates. For the three reference genes, 
geometric means (normaUzation factor: NFn, n = 3, [95]) of Ct 
and E values were calculated, and these geometric means were 
next used to calculate relative expression (R) of target genes. 
Relative gene expression levels (R) were calculated has(;d on the 
methods described in [96], in whic:h amplification (efficiencies (E) 
for target and reference samples are not assumc-d to be equal 

W) = 'E^targepT'''^'" f.»™/™/-,sr™i*)/g^^^C^,r«/ (control-sample)^ ^^^^^ 

target refers to the target gene under investigation and ref denotes 
the geometric mean of C x- and E values for the three reference 
genes. When calculating R values for constitutive expression, 
control refers to a non-template control, and sample refers to a 
biological control (i.e., untreated termites). However, when 
calculating R values for induced target gene expression, control 
refers to untreated controls, and sample refers to a treatment type 
(i.e., 0.5% PB or E. coli/P. termitis). 

R values were compared between constitutively expressed and 
induced target amplicons within and among 1) three colonies, 2) 
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